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Positions and sizes of X-ray solar flare sources 
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ABSTRACT 

Aims. To investigate the positions and source sizes of X-ray sources taking into account Compton backscattering (albedo). 
Methods. Using a Monte Carlo simulation of X-ray photon transport including photo-electric absorption and Compton scattering, 
we calculate the apparent source sizes and positions of X-ray sources at the solar disk for various source sizes, spectral indices and 
directivities of the primary source. 

Results. We show that the albedo effect will alter the true source positions and substantially increase the measured source sizes. The 
source positions are shifted up to ~ 0.5" radially towards the disk centre and 5 arcsecond source sizes can be two times larger even for 
an isotropic source (minimum albedo effect) at 1 Mm above the photosphere. X-ray sources therefore should have minimum observed 
sizes, thus FWHM source size (2.35 times second-moment) will be as large as ~ 7" in the 20-50 keV range for a disk-centered point 
source at a height of 1 Mm (~ 1.4") above the photosphere. The source size and position change is the largest for flatter primary 
X-ray spectra, stronger downward anisotropy, for sources closer to the solar disk centre, and between the energies of 30 and 50 keV. 
Conclusions. Albedo should be taken into account when X-ray footpoint positions, footpoint motions or source sizes from e.g. 
RHESSI or Yohkoh data are interpreted, and suggest that footpoint sources should be larger in X-rays than in optical or EUV ranges. 
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Introduction 

Hard X-ray (HXR) emission produced via collisional 
bremsstrahlung from solar flares by non-thermal electrons 
is the primary diagnostic tool of electron acceleration and 
transport. The spectral and spatial distributions of HXR sources 
provide us with vital clues to improve current understanding of 
the underlying physics involved in energetic electron accelera- 
tion and tra nsport. While recen t (Hard X-ray Telescope (HXT) 
on Yohkoh. iKosugi et al.lll99ll) . and mod ern (Reuven R amaty 
High Energy Solar Spectroscopic Imager. ILin et alJl20 Q2). solar 
HXR telescopes have provided superb X-ray image resolution, 
indirect imaging using either pairs of occultation grids on 
Yohkoh or rotating modulating collimators on RHESSI have (i) 
limited dynamic range and (ii) often provided inadequate spatial 
resolution exceeding the size of e.g. EUV footpoints or vertical 
extend of the chromosphere. Thus RHESSI image resolution for 
the majority of solar flares is limited to about 7", while the solar 
chromosphere is only about 3" thick. Nevertheless, unprece- 
dented spatial measurements can and have been achieved using 
the moments of X-ray distributions. The total flux (zeroth mo- 
ment) from individual sourc es in various energy r anges allows 
imaging spectroscopy (e.g. iKrucker & Linll2002t lEmslie et al] 
120031: iBattaglia & Benz) 120071: ISaint-Hilaire et alj 120081) . The 
measurements of the X-ray source positions (first moments) pin- 
point source locations with 1" or better acc uracy and allow us 
to infer the chromospheric density structu re dAschwanden et all 
l2002tlLiu et alJ200filKontar et al.l2008bh . The motions of HXR 
footpoint locati ons have been used to infer the reconnection rate 
in solar flares dFletcher & Hudson! l2002t IKrucker et al]j2003 
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Fivia n et al ]l2009ir U sing X-ray vis i bilities dHurford et all2 002 ; 



Schmahl et al. 20071) iKontar et ail (l2008bl) have measured not 
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only the positions but the HXR footpoint sizes (second moment) 
at various energies and heights and found that HXR sources 
decrease wit h energy and co nsequently with height above the 
photosphere. lXu et al.l ((2008) have measured coronal sources to 
infer acceleration region sizes. HXR images can also be inverted 
(e.g. lBrown et a l. 2006) to find the spatial electron di stributions 
and h ence the locations of electron centtoids (e.g. Prat o et al.l 
120091) . 

Since the solar atmosphere above HXR sources is optically 
thin, X-rays are often directly related to the emitting electrons. 
However, the photons emitted downwards, toward the denser 
layers of the atmosphere interact with free or bound electrons 
and can also be scattered toward the observer JTomblidll972l 
Santangelo et all fl973l) . Photons back-scattered and emerging 
back from the dense solar atmosphere to the observer create 
the albedo X-ray photons. Even for an isotropic X-ray source 
(the minimum albedo), the albedo flux can account for up to 
40% of the detected flux in the range between 30 and 50 keV 
dBai & Ramatvlll978MZhang & Huang1l2004l:lKontar et all200d 
Kasp arova et al.ll2007l) . Therefore, all X-rav sources at the so- 
lar disk are viewed as a combination of both the primary and 
backscattered fluxes. Accounting for the albedo effect is impor- 
tant for all X-ray solar observations, which can only view disk 
sources as a combination of the primary photon flux and the 
backscattered photon flux. The backscattered component taints 
the primary source properties such as electron angular, energy, 
and spatial distributions. Albedo changes the shape of the spa- 
tially integrated X-ray spectrum, which is flattened at lower en- 
ergies up to around 20-30 keV and can even produce a rtificial 
spectral features in observed spectra (Kontar et al. 2008a), while 
at higher energies above around 70 keV, the the spectrum is 
steeper than the emitted (primary) specnum. [Kontar et al.l (f2006) 
have developed and implemented albedo correction for spec- 
tral X-ray RHESSI analysis using Green's functions approxi- 
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mations by iMagdziarz & Zdziarskil (119951) . Since the reflected 
X-rays come from a rather large area (albedo patch), the surface 
brightness of the alb edo patch at the solar surface is rather low 
dBai & Ramatv fl978l) . Thi s fact explains the difficulty in directly 
imaging the albedo patch dSchmahl & Hurfo rd 2002), but high- 
lights the importance of the inclusion of albedo for understand- 
ing the measurements of the source positions and sizes (first and 
second moments), the quantities which are integrated over the 
full area of the source. 

In this Letter, using Monte Carlo simulations of X-ray 
photon transport we demonstrate how the observed positions 
and source sizes are affected by the albedo effect for various 
anisotropics, primary source sizes and primary source spec- 
tra and show that on-disk HXR sources should have energy- 
dependent minimum observed sizes. 



Spatial characteristics of the primary, backscattered 
and observed X-ray distributions 

The backscattered flux and albedo effect are studied using a 
Monte Carlo simulation starting with a hundred million pho- 
tons per run. An unpolarized X-ray source was modeled in 
space with a 2-dimensional circular Gaussian ~ exp[-jt 2 /(2ii 2 )- 
y 2 /(2d 2 )] with width d, placed at the height h — 1 Mm above 
the photosphere [The photosphere is defined here as a layer 
with H ydrogen number density 1.16 x 10 17 cm 3 (IVernazza et alJ 
fl98l l. Th is is the typical hard X-ray source height fou nd in 
footpoints dAschwanden et al J 120021: iKontar et all [2008b). The 
energy spectrum for photons has a power law 1(e) ~ e~ y with 
a spectral index of y, for energies between 3 keV and 300 keV, 
typical for RHESSI. The code accounts for the curvature of the 
Sun and the photons are assumed to move freely until they reach 



the photospheric density at a height z Q = -^R% - x 2 - y 2 - R , 

where R Q = 6.96 x 10 10 cm is the solar radius. Below this 
level photons can be either scattered or photo- electrically ab- 
sorbed. Similar to previous MC s imulations (Bai & Ra matvl 
[1971 iMagdzia rz & Zdziarskil fl995h . the Klein-Nishina cross- 
section for unpolarized X-ray radiation was used 



do ~c, 1 2 

^Q (eoA) =2 r ° 



\3 I \ 2 1 

- +- - -I sin 2 0, 



(1) 



where eo is the initial photon energy, e is the new photon en- 
ergy, 6 S is the angle between the initial and new photon direction 
and ro = 2.82 x 10~ 13 cm is the classical electron radius. After 
a scattering, the new photon energy is just given by e = eo/(l + 
-^j (1 - cos 9 S )). The absorption of X-ray photons, which is the 
dominant process below ~ 10 keV was modeled us ing modern 
solar phot ospheric abundances (|A.splund e~a l. 2009) and cross- 
sectio ns dHenke et al.lll982t iBalucinska-Church & McCammonl 
119921) for the most important elements H, He, C, N, O, Ne, 
Na, Mg, Al, Si, S, CI, Ar, Ca, Cr, Fe and Ni. For X-ray en- 
ergies > 10 keV, photoelectric absorption was approximated 



as cr fl (e () ) ~ e Q 3 (Magdzia rz & Zdziarskil 1 1995b . To account 
for elements with more electrons than Hydrogen, e.g. Helium, 
Carbon etc, Equation dT|i was multiplied b y 1.18. Our simula- 
tions differ from previous sim ulations (e.g. lBai & Ramatvll 19781: 
Magdziarz & Zdzia rskil 19951) because of newer abundances and 
the inclusion of the curvature of the Sun. The escaping photons 
are accumulated to create the brightness distribution I(x, y) over 
a given energy and solid angle. The total primary or reflected flux 
is then just an integral over the corresponding area J I(x, y)dxdy. 



Fig. QJ shows the primary and escaping photon brightness dis- 
tributions for a s ource located at the d isk centre. Similar to the 
previous results dBai & Ramatvl ll978) we see that for a com- 
pact primary source of size d — 1 .5h, the back-scattered (albedo) 
photons are reflected from an area much larger than the primary 
source. The reflected photons change the spatial distribution of 
the observed photons and produce a halo around the primary 
source. Importantly, even a primary point source will be seen as 
a source of finite size (Fig. [2]). The brightness distribution of a 
large primary source of d = 4.5h is less influenced by the re- 
flected photons but nevertheless the source will look larger than 
it actually is. 

Using solar disk centered coordinates, the centroid position 
of the source (x, y) can be found by calculating the first normal- 
ized moment of the distribution (mean) 



XCO 
oo xI(x,y)dxdy 

XCO 
co I(x,y)dxdy 



XCO 
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and the normalized variance of the distribution (second mo- 
ment), 



pOO ry pCC rt 

LJ- X ~ *) 1 (x > y) dxd y 2 J-Jy ~ y) 1 ( x ' y) dxd y 
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(3) 



Hereaf t er, following RHESS I measurements dKontar et al.l 
I2008bl: iDennis & Pernakl 120091: iPrato et al.l 120091) we will use 
source sizes in terms of FWHM (Full Width Half Maximum), 
FWHM x ,y = 2 V21n2o-.„.. 

The scattered X-ray flux depends on the cosine of heliocen- 
tric angle of the source (p = cos(0)) or equivalently on the po- 
sition of the source at the solar disk, /i = J I — (x 2 + y^/R 2 ^. A 

circular X-ray source located above the centre of the disk will 
produce a circular albedo patch (Fig.QJ). Naturally, the location 
of the HXR source and albedo patch will coincide at the disk 
centre, so albedo will not change the source position. However, 
the albedo will make the source larger than it is actually is (Fig. 
QJ). The albedo contribution becomes asymmetric if the source 
is located away from the disk centre at a given heliocentric angle 
6 (Fig. QJ-d). Due to the spherical symmetry of the Sun, there 
are two distinct directions: radial - along the line connecting 
the centre of the Sun and the X-ray source r, and perpendicu- 
lar to the radial r ± . There is no change in centroid position in the 
redirection for a spherically symmetric primary source. In the 
r-direction, the albedo causes a centroid shift towards the disk 
centre that rises from ~ at fj. — 1.0 and peaks shortly before 
falling to ~ again at fi — 0.0. Fig.[T]also shows how the source 
size varies in r ± direction, with the FWHM of the source gen- 
erally decreasing at lower /i. In the radial direction, the FWHM 
of the total and primary sources decreases close to linear due 
to a simple projection effect. The detailed 3D structure of the 
source is required before any physically meaningful predictions 
can be made concerning the change in source size in the radial 
direction, and this is beyond the scope of the paper. Therefore, 
we consider the source sizes in the r ± direction and the source 
position in the radial direction rather than along the East- West 
and So uth-North direction s. Similar to the spatially integrated 
albedo dKontar et al.ll2.006h . the shift in centroid position and the 
growth of the source are also energy and /i dependent. In the 
following, we consider the position and source size changes for 
various a) spectra of the primary source, b) primary source size, 
and c) X-ray directivity (the ratio of downward to upward emit- 
ted photons) separately. The results are summarized in Fig. [2] 
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Fig. 1. The X-ray scatter distributions of the primary photons (red dots) and the Compton back-scattered photons (blue dots) for a 
primary source at h — 1.0 Mm with d = 1.5 Mm (FWHM~ 4.9") between 20 and 50 keV for four viewing angles given by p. The 
yellow and green ellipses show the FWHM sizes for the primary and combined sources respectively. 
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Fig. 2. Spectral index dependency (panels a-d): The source position shift is in the radial direction due to albedo (source p = 0.55) and 
source size FWHM is in the perpendicular to radial direction for various spectral indices y for an isotropic source with FWHM~ 
4.9": green- y — 4; red - y = 3; purple -y = 2. Primary source size dependency (panels e-h): Isotropic primary source with 
y = 3: orange - point source, red - FWHM~ 4.9", blue - FWHM~ 14.6". Anisotropy dependency (ratio of downward to upward 
directed fluxes) (panels i-1): Simulations are for the primary source with FWHM~ 4.9" and spectral index y — 3: red - anisotropy=l 
(isotropic), blue - anisotropy=2, green - anisotropy=5. Graphs as a function of energy are for^u = 0.55 and the graphs as a function 
of p are for energies between 20 and 50 keV. 



Spectral index (Fig. |2^-d) - Similar to the spectral results, 
the albedo contribution from a smaller spectral index produces 
the largest shift in position and the larger total source size (Fig. 
|2^-d). An isotropic source of FWHM~ 4.9" with the smallest 
modeled spectral index of y — 2 produces the greatest shift of 
~ 0.5" at p = 0.5 - 0.6 and ~ 30 keV. This spectral index 
also produces the largest source size and has a FWHM~ 9.5" 
at p = 1 .0, compared with the other spectral indices of y = 3,4 
modeled. 

Primary source size (Fig.|2^-h) - For a fixed spectral index 
of y — 3, all primary source sizes produce the same shift in 
centroid position. The maximum shift in position occurs at p = 
0.5 - 0.6 and ~ 30 keV for all sources (Fig.|2^,g). Although the 



FWHM of the total source grows with increasing primary size, it 
is observed that the relative size of the total to the primary source 
is smaller for a larger primary source. This indicates that a larger 
primary source should have a smaller relative size increase due 
to albedo. Even an initial point source produces a total source 
with a FWHM peaking around 7" (Fig.|2]F,h). 

Anisotropy (Fig. |2j-l) - The shift in centroid position is 
larger for a higher initial downward anisotropy (the ratio of 
downward flux to upward flux) for all p and energies (Fig.|2},k). 
All shifts follow the general trend and tend towards zero at the 
centre (p = 1.0) and the limb (p = 0.0). Using y = 3 and a pri- 
mary source of FWHM~ 4.9", a directivity of 5 produces a peak 
difference of ~ 0.9" and even an isotropic source produces a 
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peak difference of ~ 0.4". The shift in source position peaks near 
/i = 0.4 - 0.6 and ~ 30 keV for a downward anisotropy of 2 and 
an isotropic source, but the shift peaks at a lower /j = 0.4 - 0.5 
for a downward directivity of 5. The stronger downward beam- 
ing of the primary source also leads to larger apparent source 
sizes for all /i and energies (Fig.[2],l). It should be observed that 
the total FWHM produced for a directivity of 5 peaks at fi ~ 0.15 
(Fig.[2p) giving an apparent FWHM~ 13". Since the fraction of 
reflected photons reduces with fi the FWHM in perpendicular 
direction can be expected to slowly decrease from disk centre 
to limb, but the FWHM actually increases, peaks at fi ~ 0.15 
and only then starts to decrease. This effect is due to the angu- 
lar dependence of the Compton cross-section. The cross-section 
is anisotropic and peaks at 90°, which allows larger number of 
photons to scatter into an observer direction for near the limb 
flare. It is this anisotropy in the scattering of the photons that 
causes the FWHM to peak at an angle smaller than fx=l.0. The 
high photon flux from a downward directivity of 5 allows the 
observation of this effect most clearly (Fig.|2j). 

Discussion and conclusions 

The results of the simulations show that albedo can substantially 
affect the precise position and source size measurements of X- 
ray sources. Therefore, the effect of albedo should always be 
(probably with the exception of limb/occulted flares) considered 
when the sizes or positions of X-ray sources are analyzed. The 
albedo displacement of the source position is radially directed 
towards the disk centre and depends on anisotropy of X-ray ra- 
diation, the X-ray source size and the spectral index of the pri- 
mary source. Similar to total reflected flux, the displacement of 
HXR source position is energy dependent. The largest displace- 
ment can be observed in the range between 30-50 keV at fi ~ 0.5 
(heliocentric angle ~ 60"). The shift in centroid position in this 
energy range is 0.1 - 0.5" for an isotropic (minimum albedo) 
source 1 .4" above the photosphere and this can be up ~ 0.9" for 
downward beaming with factor of 5. Because of the albedo, X- 
ray source sizes will be energy dependent, larger in the perpen- 
dicular to radial direction, and elliptical even for a spherically 
symmetric primary source. In the perpendicular to radial direc- 
tion, the largest growth in source size appears for sources close 
to the solar disk centre, in the energy range between 30-50 keV, 
where albedo is the strongest. Thus, an isotropic primary source 
with FWHM ~ 4.9" at 1.4" above the photosphere will have an 
apparent FWHM size of ~ 9" in the energy range 20-50 keV 
for sources in the wide range of heliocentric angles from 0° to 
~ 80°. 

The simulations demonstrate that X-ray sources will have a 
minimum size. An isotropic point source at 1.0 Mm above the 
photosphere will be measured by RHESSI as a source with a 
FWHM size of ~ 7" across. This result can expl ain larger X-ray 
footp o int sizes than EUV or optical ones (e.g. Kasp arova et al.l 
120051) . Id ennis & Pernakl d2009l) reported that the average semi- 
minor axis of 18 double source flares is about 4", while a few 
of the X-ray source sizes were found to be consistent with line 
sources along the flare ribbons. While the quantitative compari- 
son with the RHESSI observations requires additional work, we 
note that zero sizes are either the artifacts of the algorithms used 
or are due to very low source heights. 

The energy dependent character of albedo predicts that the 
source size as measured by RHESSI should grow with energy 
from 10 keV up to ~ 30 keV. Considering a large primary source 
of 14.6" across, e.g. a flaring loop, we fin d that the sourc e will 
grow up to ~ 18" at ~ 30 keV. Noteworthv. lXu et al.l(l2008l) have 



found that coronal source sizes are growing with energy along 
both the field lines and across. While the field line increase along 
the lines could be an indicator of electron transport or of the 
acceleration region size, the cross-field increase remains unex- 
plained, but is consistent with the growth of the source size due 
to the albedo. We note that the spatial changes of X-ray sources 
due to albedo have a great diagnostic potential for purely known 
anisotropy of energetic electrons. 
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